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Within constrained minimal-flavor-violation the large destructive flavor-changing Z-penguin man- 
aged to survive eradication so far. We give a incisive description of how to kill it using the precision 
measurements of the Z — > bb pseudo observables. The derived stringent range for the non-standard 
contribution to the universal Inami-Lim function C leads to tight two-sided limits for the branching 
ratios of all Z-penguin dominated flavor-changing K- and B-decays. 
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I. INTRODUCTION 

The effects of new heavy particles appearing in exten- 
sions of the standard model (SM) can be accounted for 
at low energies in terms of effective operators. The un- 
precedented accuracy reached by the electroweak (EW) 
precision measurements performed at the high-energy 
colliders LEP and SLC impose stringent constraints on 
the coefficients of the operators entering the EW sector. 
Other severe constraints came in recent years from the 
BaBar, Belle, CDF, and D0 experiments and concern 
extra sources of flavor and CP violation that represent a 
generic problem in many scenarios of new physics (NP). 
The most pessimistic but experimentally well supported 
solution to the flavor puzzle is to assume that all fla- 
vor and CP violation is governed by the known structure 
of the SM Yukawa interactions. In these minimal-flavor- 
violating (MFV) [1 , 2, 3 models correlations between cer- 
tain flavor diagonal high-energy and flavor off-diagonal 
low-energy observables exist since, by construction, NP 
couples dominantly to the third generation. In order to 
simplify matters, we restrict ourselves in the following to 
the class of constrained MFV (CMFV) [4] models, i.e., 
scenarios that involve only SM operators, and thus con- 
sider just left-handed currents. 



II. GENERAL CONSIDERATIONS 

That new interactions unique to the third genera- 
tion can lead to an intimate relation between the non- 
universal Zb L b L and the flavor non-diagonal Zd J L d z L ver- 
tices has been shown recently in [5] . Whereas the former 
structure is probed by the ratio of the Z-boson decay 
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width into bottom quarks and the total hadronic width, 
R®, the bottom quark asymmetry parameter, Ab, and the 
forward-backward asymmetry for bottom quarks, Ap^, 
the latter ones appear in many K- and B-decays. 

In the effective field theory framework of MFV [3] , one 
can easily see how the Zb L b L and Zd 3 L d\ operators are 
linked together. The only relevant dimension-six contri- 
butions compatible with the flavor group of MFV stem 
from the SU(2) x U(l) invariant operators 

1 =i(Q L Y u Y^ fl Q L )tfD»<p, 

_ (!) 

that are built out of the quark doublets Ql, the Higgs 
field cj>, the up-type Yukawa matrices Yjj, and the SU(2) 
generators r a . After EW symmetry breaking, O1.2 are 
responsible for both the effective Zb L b L and Zd{d l L vertex. 
Since all up-type quark Yukawa couplings except the one 
of the top, y t , are small, one has {YjjY^ji ~ y 2 V t *Vu 
and only this contribution matters in Eq. (nl . 

Within the SM the Feynman diagrams responsible for 
the enhanced top correction to the Zb L b L coupling also 
generate the Zd J L d l L operators. In fact, in the limit of in- 
finite top quark mass the corresponding amplitudes are 
up to Cabibbo-Kobayashi-Maskawa (CKM) factors iden- 
tical. Yet there is a important difference between them. 
While for the physical Z — > bb decay the diagrams are 
evaluated on-shell, in the case of the low-energy Z — » d? d l 
transitions the amplitudes are Taylor-expanded up to ze- 
roth order in the external momenta. As far as the mo- 
mentum of the Z-boson is concerned the two cases cor- 
respond to q 2 = Mj and q 2 = 0. 

The general features of the small momentum expansion 
of the one-loop Z -— > bb vertex can be nicely illustrated 
with the following simple but educated example. Con- 
sider the scalar integral 

Co = —2 IjTlTrT' A = (l+w) a (2) 

in 2 J D1D2-D3 
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with p 3 = 0. Note that we have set the space-time di- 
mension to four since the integral is finite and assumed 
without loss of generality 7713 7^ 0. 

In the limit of vanishing bottom quark mass one has 
for the corresponding momenta p\ = p\ = 0. The small 
momentum expansion of the scalar integral Co then takes 
the form 

00 / q 2 \ " 

n=0 v A/ 

with q 2 — (pi — P2) 2 — —2p\-p 2 . The expansion coeffi- 
cients a n are given by [B] 

(-1)" ^ fn\ x\ d l d n 

where 

1 /xilnxi x 2 lnx 2 \ 

g(xi,x 2 ) = , (5) 

x\ - x 2 \ 1 - xi 1 — X2 J 

and Xi = m 2 /m 2 . Notice that in order to properly gen- 
erate the expansion coefficients a n one has to keep Xi 
and x 2 different even in the zero or equal mass case. The 
corresponding limits can only be taken at the end. 

To illustrate the convergence behavior of the small mo- 
mentum expansion of the scalar integral in Eq. ^ for 
on-shell kinematics, we confine ourselves to the simpli- 
fied case mi — m 2 = M and 7713 = m t . We define 




for n — 1,2,.... The M-dependence of the relative de- 
viations 5 n is displayed in Fig.[T] We see that while for 
M < 50 GeV higher order terms in the small momen- 
tum expansion have to be included in order to approxi- 
mate the exact on-shell result accurately, in the case of 
M > 150 GeV the first correction is small and higher 



order terms arc negligible. For the two reference scales 
M = {80, 250} GeV one finds for the first three relative 
deviations 5 n numerically +9.3%, +1.4%, and +0.3%, 
and +1.1%, +0.02%, +0.00004%, respectively. 

Of course the two reference points M — {80, 250} GeV 
have been picked for a reason. While the former describes 
the situation in the SM, i.e., the exchange of two pseudo 
Goldstone bosons and a top quark, the latter presents a 
possible NP contribution involving besides the top, two 
heavy scalars. The above example indicates that the dif- 
ferences between the Zb L b L form factor evaluated on-shell 
and at zero external momenta are in general much less 
pronounced in models with new heavy degrees of free- 
dom than in the SM. Given that this difference amounts 
to around —30% in the SM [7J, it is suggestive to as- 
sume that the scaling of NP contributions to the non- 
universal Zb L b L vertex is in general under ±10%. This 
model-independent conclusion is well supported by the 
results of the calculations of the one-loop Zb L b L vertices 
in popular CMFV models presented in [5]. 



III. MODEL CALCULATIONS 

The above considerations can be corroborated in an- 
other, yet model-dependent way by calculating explicitly 
the difference between the value of the Zd J L d l L vertex form 
factor evaluated on-shell and at zero external momenta. 
In [5] this has been done in four of the most popular, 
consistent, and phenomenologically viable scenarios of 
CMFV, i.e., the two-Higgs-doublet model (THDM) type 
I and II, the minimal-supersymmetric SM (MSSM) with 
MF V [2] , all for small tan /3, the minimal universal ex- 
tra dimension (rnUED) model [S], and the littlest Higgs 
model [H] with T-parity (LHT) [TD] and degenerate mir- 
ror fermions [IT]. Examples of diagrams that contribute 
to the Z — > d? d % transition in these models can be seen 
in Fig.[2j In the following we will briefly summarize the 
most important findings of |5J. 
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FIG. 3: Relative difference 5Cnp in the THDMs, the MSSM, the 
mUED, and the LHT model as a function of M%, Mi- , 1/R, and 



Regions in the M? —5C-\ 



plane where |ReCMSSM(<7 



x L . 

0)| amounts to at least 2%, 4%, and 6% of |ReCsM(q 2 = 0)| are 
indicated by the red (gray), green (light gray), and blue (dark 
gray) points, respectively. In the case of the LHT model the 
shown curves correspond, from bottom to top, to the values 
/ = 1, 1.5, and 2TeV of the symmetry breaking scale. See text 
for details. 



In the limit of vanishing bottom quark mass, possible 
non-universal NP contributions to the renormalized off- 
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FIG. 4: Constraints on AC? 9 and AC within CMFV that follow 
from a combination of the Z — » bb POs with the measurements 
of B — > X s 7 and B — > X a l + l~ . The colors encode the frequen- 
tist 1 — CL level and the corresponding 68% and 95% probability 
regions as indicated by the bars on the right side of the panels. 
See text for details. 



shell Zd? L d l L vertex can be written as 



ji 



Gpe_ 
V2 71-2 



M^V^VuCMQ 2 )^^^ , (7) 

i 7^ j in the flavor diagonal and 
e, s w , and c w denote the Fermi 



where i = j = b and 
off-diagonal case. Gf, 
constant, the electromagnetic coupling constant, the sine 
and cosine of the weak mixing angle, respectively, while 
Vij are the corresponding CKM matrix elements. 

As a measure of the relative difference between the 
complex valued form factor C^p(q 2 ) evaluated on-shell 
and at zero momentum we introduce 



<5C N p = 1 



ReC NP (Q 2 =0) 



(8) 



ReC N p(g 2 = M!) 

The dependence of <5Cnp on the charged Higgs mass 
, the lighter chargino mass , the compactification 
scale 1/R, and Xl which parameterizes the mass of the 
heavy top T + is illustrated in Fig. [3] The allowed param- 
eter regions after applying experimental and theoretical 
constraints are indicated by the colored (grayish) bands 
and points. 

In the THDMs, the mUED, and the CMFV version 
of the LHT model the maximal allowed suppressions 
of ReCNp(<? 2 = Mj) with respect to ReCNp(g 2 = 0) 
amounts to less than 2%, 5%, and 4%, respectively. This 
feature confirms the general argument presented in the 
last section. The situation is less favorable in the case 
of the CMFV MSSM, since 5C M ssm frequently turns out 
to be larger than one would expected on the basis of 
the model-independent considerations if the masses of 
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TABLE I: Bounds for various rare decays in CMFV models at 95% probability, the corresponding values in the SM at 68% and 
95% CL, and the available experimental information. See text for details. 





Observable 


CM1 V (95% CL) 


bM (68% CL) 


bM (95% CL) 


Experiment 




B(K + 


— > ir + vv) x Ifl 11 


[4 29 19 721 


7 15 + 1 28 


[5 40 9 111 


(14 7+13-01 RTl 
(J^- ' _8.9 j l^ij 




B(K L 


-> Tr'W) X 10" 


[1.55,4.38] 


2.79 ±0.31 


[2.21,3.45] 


< 2.1 x 10 4 (90% CL) 


m 


B(K L - 


-> /i + /tt~)sD X 10 9 


[0.30, 1.22] 


0.70 ±0.11 


[0.54,0.88] 






B(B 


-» X d vv) x 10 6 


[0.77, 2.00] 


1.34 ±0.05 


[1.24,1.45] 






B(B 


-> X 3 vv) x 10 5 


[1.88,4.86] 


3.27 ±0.11 


[3.06,3.48] 


< 64 (90% cl) rja 




B(Bd 


-> fi+fi') x 10 10 


[0.36,2.03] 


1.06 ±0.16 


[0.87, 1.27] 


< 3.0 x 10 2 (95% CL) 


M 


B(B S 


-> a* + aO x 10 9 


[1.17,6.67] 


3.51 ±0.50 


[2.92,4.13] 


< 9.3 x 10 1 (95% CL) 





the lighter chargino and stop both lie in the hundred 
GeV range. However, the large deviation <5Cmssm & re 
ultimately no cause of concern, because |ReCMSSM(<Z 2 = 
0)/ReCsivi(9 2 = 0)| itself is always below 10%. In conse- 
quence, the model-independent bounds on the NP con- 
tribution to the universal Z-penguin function that will 



be derived in the next section do hold in the case of 
the CMFV MSSM. More details on the phenomenologi- 
cal analysis of <5C N p in the THDMs, the CMFV MSSM, 
the mUED, and the LHT model including the analytic 
expressions for the form factors Cnp(? 2 ) can be found in 
the recent article [5]. 



IV. NUMERICAL ANALYSIS 

Using the technique of cpsilon parameters a model- 
independent numerical analysis of AC = Re C(q 2 = 
0) — ReCsM(<Z 2 = 0) is a back-on-the-envelope calcu- 
lation. The variation e^ p = e& — ef M arising from NP 
contributions to Zb L b L can be defined through the inclu- 
sive partial width of Z — ► bb as follows [El] 



-pNP 
1 bb 



where 



9± 
9 A 



-- (V2G F M 2 z )i (g h v {h^) - g b A {h^ 5 b))z» , (9) 
(l + ^ + ^ P )^, A = (l+$ P )9 d A - (10) 



From Eqs. Q, and (|9j one obtains 



AC 



V2G F M*cl 



(1±<5C NP )4 



NP 



(11) 



By combining experimental 13] and theoretical uncer- 
tainties |15j in €b and linearly one finds 

£ np = (0.4 ±2.5) x 1CT 3 . (12) 

Assuming <5Cnp = ±0.1 one then arrives at 



AC = -0.04 ±0.26, 



(13) 



which implies that large negative contributions that 
would reverse the sign of the SM Z-penguin ampli- 
tude are highly disfavored in CMFV scenarios due to 
the strong constraint from [5J. Interestingly, such 



a conclusion cannot be drawn by considering only flavor 
constraints [14], since a combination of B(B — > X s j), 
B(B -> X s l + l~), and B(K+ -> tt+w) does not allow to 
distinguish the SM solution AC = from the wrong-sign 
case AC » —2 at present. 

The result in Eq. (13 1 agrees amazingly well with the 
numbers of a thorough global fit to the POs R®, Ab, 
and Ap B jT3] and the measured B — > X s j [T7] and 
B — > X s l + l~ [T5] BRs obtained by employing customized 
versions of the ZFITTER [T5] and the CKMfitter pack- 
age [19]. Neglecting contributions from EW boxes these 
bounds read [5] 



AC = -0.026 ± 0.264 (68% CL) 
AC = [-0.483, 0.368] (95% CL) 



(14) 



The constraint on AC within CMFV following from the 
simultaneous use of R°, Ab, Ay B , B(B — > X s j), and 
B(B -> XJ+1-) can be seen in Fig.g 

One can also infer from this figure that two regions, 
resembling the two possible signs of the amplitude A(b — > 
sj) oc Cy ff (mh), satisfy all existing experimental bounds. 
The best fit value for ACf = Cf(m b ) - C$% M (m b ) is 
very close to the SM point residing in the origin, while 
the wrong-sign solution located on the right is highly 
disfavored, as it corresponds to a B(B — ► X s l + l~) value 
considerably higher than the measurements [20 . The 
corresponding limits are [5] 



ACf = -0.039 ± 0.043 (68% CL) , 
ACf = [-0.104,0.026] U [0.890,0.968] (95% CL). 



(15) 



Similar bounds have been presented previously in [14] . 
Notice that since the SM prediction of B(B — > X s j) 
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[IB] is now lower than the experimental world average 
by 1.2 a, extensions of the SM that predict a suppres- 
sion of the b — > 57 amplitude are strongly constrained. 
In particular, even the SM point ACf = is almost 
disfavored at 68% CL by the global fit. 

The stringent bound on the NP contribution AC given 
in Eq. ( 14 1 translates into tight two-sided limits for the 



V. CONCLUSIONS 



BRs of all Z-penguin dominated flavor-changing K- and 
B-decays as shown in Tab.|l] A strong violation of any of 
the bounds by future measurements will imply a failure 
of the CMFV assumption, signaling cither the presence 
of new effective operators and/or new flavor and CP vi- 
olation. A way to evade the given limits is the pres- 
ence of sizable corrections <5Cnp and/or box contribu- 
tions. While these possibilities cannot be fully excluded, 
general arguments and explicit calculations indicate that 
they are both difficult to realize in the CMFV framework. 



R.I. P. large destructive CMFV Z-penguin! 
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